Theory of random matrices with strong level confinement: 
orthogonal polynomial approach 
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Strongly non-Gaussian ensembles of large random matrices possessing unitary symmetry and 
logarithmic level repulsion are studied both in presence and absence of hard edge in their energy 
spectra. Employing a theory of polynomials orthogonal with respect to exponential weights we 
calculate with asymptotic accuracy the two-point kernel over all distance scale, and show that in 
the limit of large dimensions of random matrices the properly rescaled local eigenvalue correlations 
are independent of level confinement while global smoothed connected correlations depend on con- 
finement potential only through the endpoints of spectrum. We also obtain exact expressions for 
density of levels, one- and two-point Green's functions, and prove that new universal local relation- 
ship exists for suitably normalized and rescaled connected two-point Green's function. Connection 
between structure of Szego function entering strong polynomial asymptotics and mean-field equation 
is traced. 
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I. INTRODUCTION 

Statistical properties of complex physical systems can 
successfully be investigated within the framework of 
random-matrix theory (RMT) [Q. It turned out to be 
quite general and powerful phenomenological approach to 
description of various phenomena in such diverse fields as 
two-dimensional gravity quantum chaos complex 
nuclei Q] and mesoscopic physics || . 

In all the realms mentioned above the physical sys- 
tems can be described with the help of different matrix 
models whose structures depend on physical properties of 
the systems involved. In the applications of RMT to the 
complex quantum-mechanical objects the real Hamilto- 
nian is rather intricate to be handled or simply unknown. 
In such situation the integration of exact equations is 
replaced by the study of the joint distribution function 
P [H] of the matrix elements of Hamiltonian H. If there 
is not preferential basis in the space of matrix elements 
(i.e. the system in question is "as random as possible" , 
and equal weight is given to all kinds of interactions) one 
has to require P [H] d [H] to be invariant under similarity 
transformation H — ► 7?. _1 H7?. with 1Z being orthogonal, 
unitary or symplectic n x n matrix reflecting the fun- 
damental symmetry of the underlying Hamiltonian. The 
general form of P [H] compatible with invariance require- 
ment is 



P[H] = Z -1 exp{-trV [H]} 



(1) 



with arbitrary V [H] providing existence of partition 
function Z. Introducing matrix Sp that diagonalizes 
Hamiltonian H, H = S7 XSp, and carrying out the 
integration over orthogonal (/? = 1), unitary (/3 = 2) or 
symplectic (/3 = 4) group d/i (Sp) in the construction 



P [H] d [H] , one obtains the famous expression for the 
joint probability density function of the eigenvalues {x} 
of the matrix H: 



P({x}) 



exp 



(2) 



Level repulsion described by the logarithmic term is orig- 
inated from Jacobian J| 







i<3 



arising when pass- 
ing from the integration over independent elements Hij of 
the Hamiltonian H to the integration over smaller space 
of its n eigenvalues {x}. Confinement potential V(x), 
which determines (together with logarithmic law of level 
repulsion) the mean level density, contains an informa- 
tion about correlations between different matrix elements 
of a random Hamiltonian H. [Note, that parameter /3 is 
factored out from V [H] in Eq. ([j]) to fix density of levels 
in random-matrix ensembles with the same confinement 
potential but with different underlying symmetries]. 

In the matrix formulation given above the eigenval- 
ues {x} of the Hamiltonian H run from — oo to +oo. 
Formally, the same matrix model Eq. (Q) appears in so- 
called maximum entropy models constructed to describe 
transport properties of mesoscopic systems. In this case 
there is additional positivity constraint on {x}, x > 0, 
that directly follows from unitarity of scattering matrix 
^,[1 and introduces the hard edge into eigenvalue spec- 
trum. 

In the unitary case (j3 = 2), which applies to physical 
systems with broken time-reversal symmetry, the struc- 
ture of Eq. (|^) allows one to represent exactly all the 
global and local statistical characteristics of the physical 
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system, such as averaged density of levels, n-point corre- 
lation functions, level-spacing distribution function etc., 
in terms of polynomials orthogonal with respect to the 
weight function w (x) — exp {—2V (x)} on the whole real 
axis R (or on R + if there is a hard edge in eigenvalue 
spectrum). [Otherwise, when (3 = 1 or /3 = 4 more com- 
plicated sets of skew orthogonal polynomials should be 
introduced 0]. 

Analytical calculation of the corresponding set of or- 
thogonal polynomials is a non-trivial problem. However, 
if the elements Hij of the random matrix H are believed 
to be statistically independent from each other, one ob- 
tains the quadratic confinement potential V (x) ~ x 2 || 
leading to the Gaussian Invariant Ensembles of random 
matrices. In such a case there are significant mathemati- 
cal simplifications allowing to solve the matrix model Eq. 
(@) completely §. 

From the very beginning it was understood || that 
requirement of statistical independence of the matrix el- 
ements is not motivated by the first principles, and, 
therefore, several attempts were undertaken to elucidate 
an influence of a particular form of confinement potential 
on the predictions of the random matrix theory developed 
for Gaussian Ensembles. 

Two essentially different lines of inquires of this prob- 
lem can be distinguished. The first line lies in the frame- 
work of polynomial approach, while a second one consists 
of developing of a number of approximate methods. The 
mean-field approximation proposed by Dyson JhJ allows 
to calculate density of levels in random-matrix ensemble. 
This approach being combined with the method of func- 
tional derivative of Beenakker makes it possible 
to compute global (smoothed) eigenvalue correlations in 
large random matrices. Smoothed correlations can also 
be obtained by diagrammatic approach of Brezin and 
Zee [ [l3| and by invoking the linear response arguments 
and macroscopic electrostatics [H. We stress that all 
the methods mentioned above allow to study correlations 
only in long-range regime and, in this sense, they are less 
informative as compared with the method of orthogo- 
nal polynomials [Q. It is worth pointing out the super- 
symmetry formalism Jl^| , recently developed for matrix 
model Eq. @) with non-Gaussian probability distribu- 
tion function P [H] , which is exceptional in that it allows 
to investigate local eigenvalue correlations and represents 
a powerful alternative approach to the classical method 
of orthogonal polynomials. 

In the framework of polynomial approach there was 
a number of studies to go beyond the Gaussian distri- 
bution P[H]. In Rcfs. @ - @ it was found out that 
unitary random-matrix ensembles associated with clas- 
sical orthogonal polynomials exhibit Wigner-Dyson level 
statistics [for corresponding ensembles with orthogonal 
and symplectic symmetry see Ref. pEj]. Non-Gaussian 
unitary random-matrix ensembles associated with (sym- 
metric) strong confinement potentials V (x) = x 2 + 7a; 4 



and V (x) — J2n=i a n x2n were treated in Refs. and 
pot , respectively. (We note that both potentials men- 
tioned above are stronger than quadratic, and they do not 
refer to the maximum entropy models). As far as these 
works have been based on different conjectures about 
functional form of asymptotics of polynomials orthogonal 
with respect to a non-Gaussian measure, and the prob- 
lem of hard edge in eigenvalue spectrum was out of their 
scope, the polynomial approach to basic problems of the 
random matrix theory needs further and more rigorous 
study. 

The purpose of the present work is to show that the 
problem of non-Gaussian ensembles with unitary sym- 
metry does can be handled rigorously by the method of 
orthogonal polynomials. Our treatment is exact (i.e. it 
does not involve any conjectures) and based on the recent 
results obtained in the theory of polynomials orthogonal 
with respect to exponential weights on R. It applies to 
very large class of confinement potentials which is much 
richer than that considered in Refs. @J|^] and allows also 
to treat the matrix models with positive constraints on 
eigenvalue spectrum. We concentrate on the calculations 
of density of levels, one- and two-point Green's functions, 
two-point kernel and connected "density-density" corre- 
lation function over all distance scale. This allows us to 
resolve the problem of universality for local and global 
correlations of random-matrix eigenvalues and to estab- 
lish a new universal local relationship for properly nor- 
malized and rescaled connected two-point Green's func- 
tion. One of the interesting points we would like to stress 
is that the mean-field approximation, widely used in 
the theory of random matrices, naturally appears in our 
treatment without any physical speculations and turns 
out to be closely allied with structure of Szego function 
entering strong pointwisc asymptotics of orthogonal poly- 
nomials. 

The paper is organized as follows. Section II contains 
a short introduction to the theory of polynomials orthog- 
onal with respect to the Freud weights. The asymptotic 
formula for orthonormal "wave function" that we need in 
later sections is given there. In Section III we calculate 
two-point kernel and resolve the problem of universality 
of level statistics. The density of levels and one-point 
Green's function are computed in Section IV. Connec- 
tion between structure of Szego function and mean-field 
equation is established there as well. Section V is devoted 
to the calculation of the two-point connected Green's 
function; corresponding new universal local expression is 
given. Section VI contains generalizations of the results 
obtained in the preceedings Sections for a wider class 
of random matrices characterized by Erdos-type confine- 
ment potential. In Section VII we present a treatment of 
the maximum entropy models with hard edge. Finally, 
in Section VIII we discuss the results obtained. 



2 



II. FREUD-TYPE CONFINEMENT POTENTIALS 
AND CORRESPONDING ORTHOGONAL 
POLYNOMIALS 

Let us consider a class of symmetric (even) confine- 
ment potentials V (x) supported on the whole real axis 
x € (—00, +00) which are of smooth polynomial growth 
at infinity and increase there at least as \x\ 1+S (S is arbi- 
trary small positive number). More precisely, we demand 
V (x) and cPV/dx 2 be continuous in x G (0, +00), and 
dV/dx > in the same domain of variable x. We also 
assume that for some A > 1 and B > 1 the inequality 



A < 1 + x 



d 2 V/dx 



< B 



(3) 



dV/dx 

holds for x G (0, +00), and also for x positive and large 
enough 



\d 3 V/di 



< const. 



(4) 



dV/dx 

The class of potentials V (x) satisfying all the above re- 
quirements is said to be of the Freud type [|TJ . The typi- 
cal examples of the Freud potentials are (i) V (x) = \x\ c 



l°ln"(7- 



with a > 1, 



with a > 1, and (ii) V (x) = \x\ 
(3 £ R, and 7 large enough. 

Now it is convenient to introduce a set of poly- 
nomials P n (x) orthogonal with respect to the Freud 
(non-Gaussian) measure day: (x) — wjr (x) dx — 
exp(-2V(x))dx, 

P n (x) P m (x) day (x) = S nm , (5) 

for which the following basic result was obtained by D. 
S. Lubinsky & 



lim 



+ 1 



d\- 



1/2 



xRe 



z n D- 



F„ 



Wjr (o n A) = 0. 



(6) 



Here parametrization z = e and A = cos 8 was used. 

Szego function D (g; z), appeared in Eq. (^J), is of fun- 
damental importance in the whole theory of orthogonal 
polynomials p3| , and takes the form 



1 f +7V 1 
D (g; z) = exp ( — / dy— 



(7) 



The first argument of Szego function in Eq. (|6|) is 

1/2 

F n (<p) = exp(-y (a n cosLp)) |sin</j| , (8) 

and a n is the n-th Mhaskar-Rahmanov-Saff number be- 
ing the positive root of the integral equation 0] 

r All ,\ I II I \ 

" = / , = I — I ■ (9) 



2a n 

-■' Jo \f\ - A 2 V dx 



(In what follows it will be seen that a n is none other 
than band edge for eigenvalues of corresponding random- 
matrix ensemble.) 

Equation (^) may be rewritten in the different form 
passing on to the integration over x — a n X (so that 
parametrization x — a n cos 8 takes place) : 



r+a„ ( 
lim / dx < P n (x) 



2 



1/2 



z n D~ 2 ( F n : - 



xRe 



Analogously, Eq. (m reads 



2 

n = — 



Wjr (x) = 0. 



cdx dV 



7T Jo y/al - x 2 dx 



(10) 



(11) 



Since from Eq. ( |ll|) it follows that limn^oo a n =/= 0, we 
immediately conclude that expression in parentheses of 
Eq. ( |Io|) asymptotically tends to zero as n — > 00 on the 
interval of integration \x\ < a n . If one is not interested 
in remainder term, we arrive at the asymptotic formula 
for orthogonal polynomials of the Freud type: 



Pn(x) = x Re 

na n 



z n D- 2 (F n ;- 



, X £ (—On, +a n ) ■ 
(12) 

Szego function D (g; e lS ) may be represented as [Q 
D(g;e w ) = y/g~{8) exp [iT (g; 6)] , (13) 

where 



"V <" 1 ( '-^r~ ) [ m 9 (V>) - In 5 ■ 

(14) 

Making use of the representation of Eqs. (|l^) and 
@, noting that F n (-ip) = F n (ip) and T (F n ; —8) = 
— r (F n ; 0), we obtain 



D[F n ;-j = exp (-^V(a n cos( 



x |sin6»| 1/4 exp[-ir(i^;l9)] . 



(15) 



Then, Eqs. (|12|) and Ql5|) yield 

Pn («n COS 8) = 



2 



exp (V (a n cos 8)) 



sin 1 



,1/2 



{n8 + T(F 2 ;8)), (16) 



where 



3 



4tt 



dip cot 



[lnF 2 (^)-lnF 2 



= — / d<p [in F^(<p)- In (9)} ^ cot 



cot 



— / ^[lnF 2 (^)-lnF 2 (fl)] 



sm( 



Introducing the new variable of integration £ = a n cos <p 
and using parametrization x = a n cos8 (|x| < a n ), we 
get 



7n (x) - r (-F, 



l x—a n cos ( 



5?' * 



y/c%=rfh(Z)-h(x) 



with 



M0 = -2^(0 + 2 ln 



1 - 



Since for |x| < a ? , 
V 



d£ 



-a» K-sJv/a^ 3 



(18) 



(19) 



(20) 



(here 7-* stands for principal value of an integral), Eq. 
( |l8| ) can be rewritten in the form 



7n (a:) - 7T V 

Z7T 



(21) 



Then we obtain following asymptotic formula for or- 
thonormal "wave functions" ip n (x) = P n (x) exp(— V (x)) 
that we need in what follows: 
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-1/4 



n arccos — + 7„ I x j 



(22) 



We remind that Eq. is valid for \x\ < a n in the limit 
n — > oo. 



III. TWO-POINT KERNEL AND UNIVERSAL 
EIGENVALUE CORRELATIONS 

Two-point kernel allowing to calculate all the global 
and local characteristics for the random-matrix ensem- 
bles is determined as M 



K n {x, y) 



1 

ir(y-x) 



K n (x, y) 



cos ip — COS t) 

k n -l (y) V'n-l (x) - 1p n (x) 1p n -l (y) 



(17) 



y-x 



(23) 



where k n is a leading coefficient of the orthogonal poly- 
nomial P n (x) . Substitution of Eq. ( p^ ) into Eq. (j2~ 
yields in the large-n limit 



K n (x, y) = 



2 krj 



7TCL n k n 



y-x 



-1/4 



x (cos $„-i (x) cos $„ (y) - cos $„_i (y) cos $„ (x)) , 

(24) 



where 



3?n (a?) = 7n (#) + n arccos 



(25) 



In Eq. ( p4| ) the fact was used that lim, woo (a n -i/a n ) — 
1. Really, as was noted in previous section, the Freud- 
type potentials exhibit a polynomial growth at infinity. 
Supposing that at large positive x potential V (x) roughly 
behaves as x p (p > 1) we immediately obtain the esti- 
mate (see Eq. (p])) a n — > n x / p as n — > oo. Then, obvi- 
ously, lim„^oo (a„_i/a„) = 1. Taking into account this 
limit and carrying out the changing of integration vari- 
able £' = £a„/a„_i in Eq. ( pl| ) we easily obtain that in 
the large-n limit j n —i (x) = j n (x), and as a consequence 



(26) 



<fr„_i (x) = (x) - arccos — 



Now Eqs. (||) and (§6) give us 



-1/4 
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x } cos$„ (z)cos$„ (y) 



sin <& n (y) cos$„ (x) 



sin$„ (x) cos$ n (y) \ 1 



(27) 



When deriving we have used identity linin^oo k n ^i/k n a n = 
1/2 proved in Ref. |||. We stress that Eq. @ is valid 
for arbitrary x and y lying within the band (— a n , +a n ). 

Equation ( p7f ) allows to determine smoothed (over the 
rapid oscillations) connected correlations v c (x, y) of the 
density of eigenvalues v n (x) p2],p0[, 



(x, y) = (v n (x) v n (y)) 



- {v n (x)) {v n (y)) = -Kl (x, y) 1 x^y (28) 

by averaging over intervals |Ax| <C a n and \Ay\ <C a n 
but still containing many eigenlevels. Direct calculations 
yield simple universal relationship 

v c (x, y) = 



K - X V 



F , x £ y (29) 



2tt (x-y) y/al - x 2 ^a 2 n - y 

with dependence on the potential V (x) only through the 
endpoint a n of the spectrum. 

Now we turn to the local properties of two-point ker- 
nel. Assuming that in Eq. ( p7| ) \x — y \ <C a n and both x 
and y stay away from the (soft) band edge a n we obtain 



K n (x, y) 



sin ($» (x) - g» (y)) 
7T (y - x) 



(30) 



where (x) is defined by Eq. (|25|). This two-point ker- 
nel may be rewritten in locally universal form. Taking 
into account the integral representation 

$n 0) = \ arccos ( — ) - w f co an (f ) d£ + ^ (2n - 1) , 
2 \a n / Jo 4 

(31) 



w «» (&) = -o^ / 75 I?T7T rs 7 *> ( 32 ) 



proved in Appendix, we see that Eq. ( |30| ) may be rewrit- 
ten as 



K n (x, y) 



sm 



7T (y - x) 



(33) 



The characteristic scale of the changing of w an (0 is 
[oj~^ du> ' an I 'dQ ~ a„, so that for |x— y| <C a„ (that 
has been supposed in Eq. (|3(])) Eq. ( |33| ) is reduced to 
universal form 



#n (as, y) 



sin \nv n (y - x)] 
tt (y - x) 



(34) 



with 77„ = Lu an (^2 ) P^ a y m S ^ ne r °l e °f l° ca l density of 
levels. Correspondingly, local two-level cluster function 
being rewritten in rescaled variables s and s' 



Y 2 (s, a') 



K 2 n (x,y) 



\v n (x)) (u n (y)) 



y=y(s') 



sin 2 [tt (s — s' 



(35) 



proves universal Wigner-Dyson level statistics in the uni- 
tary random-matrix ensemble with Freud-type confine- 
ment potentials (here s = V n x and s' = v n y are the 
eigenvalues measured in the local mean level-spacing). 



IV. DENSITY OF LEVELS AND ONE-POINT 
GREEN'S FUNCTION 

Expression for density of levels defined as 

(v n {x)) = <tr<5(x-H)) =K n (x,x) (36) 
immediately follows from Eq. pO]): 

K (x)> = — = - - - — , (37) 

tt dx tt y y/al - x 2 ax J 

(see Eq. ©). Using Eqs. (|l|) and @, and 
parametrization x = a n cos0, we obtain the formula 

1 

(z/„ (x = a n cos 0)) 



7ra„ sm t 



cZ0 



argD r 



e 2y ( a " cos ^) \smip\;e v 



nO 



(38) 



which establishes the connection between density of lev- 
els in random-matrix ensemble with Freud-type confine- 
ment potential and Szego function for corresponding set 
of orthogonal polynomials, Eq. |7|). 
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Another representation of level density can be obtained 
from Eqs. and @: 



' J n 0*0) = Z2 V 



(39) 



This formula is rather interesting and deserves more 
attention. Considering this expression as an equation for 
dV/dx one can resolve it invoking the theory of integral 
equations with Cauchy kernel pq|: 



^{X 1 )) ^ = dV 

x — x 1 dx 



(40) 



Thus, one can think that density of levels is a solution of 
integral equation 



V(x) 



+a„ 



dx' (v n (x 1 )) In \x — x | + fi (41) 



with /i being "chemical potential" . It is non more than 
famous mean- field equation which, in our treatment, fi- 
nally follows from the asymptotic formula Eq. ( |l2] ) for 
orthogonal polynomials. Quite surprisingly, Szego func- 
tion Eq. (0) turns out to be closely related to the mean- 
field approximation by Dyson |Io[| . 

Now we can easily calculate the one-point Green's 
function 



G p (x) 



tr- 



1 



x - H + ipO 



1 



(tr£(£-H)) 



' x — £ + ipO 

where p = ±1. Last integral can be rewritten as 

(Vn (0) 



(42) 



G p (x) = V 



x-i 



- i^P {v n (x)} , (43) 



whence we obtain by means of Eqs. (p9) and (Eoft 



GP(x) = dV_2ip p 
dx IT 



idi dVy/aJ^ 



(44) 



We would like to stress that both Eqs. (|39j) and ([44]) 
have been obtained within the framework of the theory of 
polynomials orthogonal with respect to the Freud mea- 
sure. This comment equally pertains to the mean-field 
equation Eq. (Ell). 



H 



tr- 



tr- 



x v — H / \ x'i — H / 



(45) 



where x p = x + ipO and x' p , = x' + ip'O (p,p' = ±1). It 
can be rewritten in an integral form 



Gf (x,x') = 



d£dr/ 



x [{v n (£) v n {rf)) - {v n (£)) {v n (v))} 



(46) 



Recognizing that the quantity in parentheses is 

{v n (0 »n {v)) c = {"n (£)) $ (£ - v) ~ K n (f , v), we obtain 
the formula 

Gf (x,x') = 



+a„ 



^ K (£)> 



TT 2 pp'K^ (x, x') 



+ iir [pA (x, x ) + p A (x , x)] — A (x, x ) , 
where the following notations were used: 

-KM 



A(x,x')=V 



x' -£ 



A (,-.,') p j +aU d£±^-. 



(47) 



(48) 



(49) 



Two-point kernel K n (x,x') entering Eqs. J47|), ( EM ) is 
determined by Eq. (p7|). 



A. Smoothed connected two-point Green's function 

Let us consider the first integral Eq. Substituting 
Eq. ( p7|) into Eq. (EJh and taking into account that terms 
of the type sin$„ (£), cos$ n (£), and sin$„ (£) cos$„ (£) 
oscillate rapidly and, therefore, do not contribute into in- 
tegral over £ in the leading order in ii>l,we have after 
some rearrangements 



A (x, x') = 



1 



1 



2^ 2 y/1 - x */al 



V. TWO-POINT CONNECTED GREEN'S 
FUNCTION 

Two-point connected Green's function is defined as 



x V 



- a „ (x'-0 (x-tf y/l-?/a* 



xZ 



(50) 



G 



provided x ^ x'. Formally, this integral is divergent 
thanks to the double pole of integrand oc (x — £)~ 2 . It 
is easy to see that this singularity is rather artificial and 
connected with the fact that condition x ^ £ was sup- 
posed to be fulfilled when neglecting rapid oscillations in 
£ in integrand of Eq. (^). This is the reason why the 
integrand in Eq. ( |50| ) displays a non-correct behavior in 
the vicinity x = £. Actually, as can be verified, the in- 
tegrand is finite for x = £, and corresponding integral 
is convergent. Moreover, direct comparison of Eq. ( |50| ) 
with results |l2] ] shows that equation in question can be 
rewritten in the form 



A(x,x')=V 



+a„ 



where 



(51) 



we finally arrive at the expression for two-point connected 
Green's function: 

1 



G pp ' (x,x>) 



x \PP 



XX 



(x - X 1 ) V a n ~ X 2 



(x p x' p , ^ 



(57) 



Here we have used Eqs. (|2S|) and (|29|). Equation (|57| ) 
is valid for arbitrary x ^ x' lying within the band 
(— a n , +a„). Universal relationships of this type were ob- 
tained in Ref. [El. 



T 2 (£, x) = K 2 (£, x) + (v n (0) S (£ - x) 
is two-level cluster function, and 

i s K (0) 



K 2 (£,x) 



2 SV(x) 



(52) 



(53) 



is two-point correlation function (the notations of Ref. 
p2| have been used). Then, taking into account Eqs. 
(pi]), ( |52| ) and (p)|), we obtain from Eq. ( |47| ) after some 
transformations: 



G pp (x, x') = ir 2 pp'Kl (x, x' ) + in 



pV 



+a „ 



d( 



x'-i 



K 2 (^x) +p'V 



VV 



d£dr/ 



K 2 (Z,r)). (54) 



Now we only have to calculate the integrals containing 
K 2 . The most proper way is to invoke the integral equa- 
tion 



and definition Eq. (|53|). Since Eqs. (|5 
tity 



(55) 

|) yield iden- 



pV 



K 2 (Z,x)+p'V 



#2 



B. Local connected two-point Green's function 

In the local regime, when \x — x'\ <C a„, one cannot 
disregard oscillations of integrands in Eqs. (|4^) and (|4S|). 
Since in this energy scale the density of states (y n {x)) is 
slowly varying function and two-point kernel K n (x, x') 
is universal, Eq. (p4|), one obtains that 28 



A (x, x' ) 



1 



sin [27ri/„ (x' — x 

(x' — x) 



and 



A (x, x') = 



sin 2 [irv n (x — x') 



(x — x') 

Then Eqs. (||), © and §T§ yield 

Gf ' (a, x') = tt 2 77„ | p - p'\ S(x- x') 



(58) 



(59) 



W - 1] 



sin [7rz/ n (x — x')] 
(x — x') 2 



, , sin [nvn (x - x)] cos [irv n (x - x)] 

+ np -p) — — ^ ■ ( 60 ) 

[x' — X) 

This equation only depends on the local mean-level spac- 
ing V n , and therefore it can be written down in univer- 
sal form. Introducing normalized and rescaled two-point 
connected Green's function 



-W 



d^drj 



a n {x ~ (x' - rf) 



K 2 (H,V) 



(x p x' p , 



(s,s') = 



G pp (x,x') 
(v n (x)) {v n (x')) 



(61) 



/(.') 



where s = v n x and s' = v n x' are the eigenvalues mea- 

(56) 

sured in the local mean level-spacing, we obtain new uni- 
versal relationship in the random-matrix theory: 
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gf (s,s')=n 2 \ p -p'\S(s-s') 

+ i(p-p') Sm kil^pi e ™(s-s'>Mp-p'). (62) 
(s - s>) 

Note that expression of this type was previously obtained 
in Ref. p9fl only for Gaussian random-matrix ensemble 
using supersymmetry formalism. 



VI. EXTENSION FOR ERDOS-TYPE 
CONFINEMENT POTENTIALS 



constraints on eigenvalues {x} (maximum entropy mod- 
els). The constraint x > is essential feature of those 
models that follows directly from the unitarity of scat- 
tering matrix and imposes the presence of the hard edge 
in the energy spectrum of the matrix model. To our 
knowledge there is no rigorous treatment of such matrix 
model with strong enough confinement potential V (x) 
within the method of orthogonal polynomials except for 
generalized Laguerre ensembles of random matrices [ |30| . 

Below we show how the problems associated with max- 
imum entropy model can be treated within the polyno- 
mial approach in very general case. 



All the results obtained above are valid for confine- 
ment potentials exhibiting smooth polynomial growth at 
infinity (see Section II) but they can be extended for an 
Erdos-type confinement potentials which grow faster than 
any polynomial at infinity (see Ref. |27|, Ch. 2). 

Namely, let V (x) be even and continuous in x S 
( — oo, +oo), d 2 V/dx 2 be continuous in x e (0, +oo), 
dV/dx be positive in the same domain of x and continu- 
ous at x = 0. Moreover, let 



T(x) 



d 2 V/dx 2 
'' dV/dx 



(63) 



be positive and increasing in x G (0, +oo) with 
linij^+o T (x) > while lim^—,,^ T (x) — oo, and 



T 



(x) = O {[dV/dx) 1,lb ^ for x 



■ oo, 



(64) 



d 2 V/dx 2 dV/dx 
dV/dx ~ V(x) 



and 



\d 3 V/dx 3 
dV/dx 



< const 



fdV/dx 



(65) 



for x large enough. The class of potentials V (x) sat- 
isfying all the above requirements is said to be of the 
Erdos-type. The simple examples of Erdos-type confine- 
ment potentials are (i) V (x) — exp fc (|a;| a ) with a > 
and k > 1 (here exp fe denotes the exponent iterated k- 
times); (ii) V (x) — exp (ln Q (7 + x 2 )) with a > 1, and 7 
large enough. 

Polynomials orthogonal with respect to the Erdos mea- 
sure das = W£ (x) dx — exp (-2V (x)) dx (here V is of 
Erdos type) have the same asymptotics p^j and, there- 
fore, Eq. (|2^) remains valid along with all the results 
obtained in Sections III, IV and V. 



VII. MATRIX MODELS WITH POSITIVITY 
CONSTRAINTS ON EIGENVALUES 



A. Polynomials orthogonal on x > 



Let confinement potential V (x) be of the Freud- or 
Erdos-type defined on the whole real axis R, that is V is 
monotonous function behaving at least as (S > 0) 

and growing as or even faster than any polynomial at in- 
finity, and P n (x) be a set of polynomials orthogonal on R 
with respect to the measure da{x) — exp {— 2V (x)} dx 
(see Eq. (^)). Then polynomials 



S n (x) = P 2 „ (y/x) 



(66) 



form a set of polynomials orthogonal on R + with the 
measure [BlJ da s (x) = exp{— 2V S (x)}dx, 



S n (x) S m (x) da s (x) = S nm , 



where confinement potential 



V 8 (x) =V (Vx) + -Inx 



(67) 



(68) 



is a monotonous function that behaves at least as |x| 3+<5 
(6 > 0) and can grow even faster than any polynomial at 
infinity. 

Equation ([36]) allows to write down large-n asymp- 
totics for introduced set of orthogonal polynomials. It 
is straightforward to get from the results outlined in Sec- 
tion II and Appendix the following asymptotic formula 
[which is analogue of Eq. jl^)]: 



S n (x) 



2exp(V s (x)) 



1 



* (x6„) 1/4 [l-x/6„] 1/4 



cos$„ (x) , 



(69) 



where x £ (0, b n ) , and 



(x) = \ arccos (J± 



In the random-matrix theory of quantum transport 
|p] the matrix model Eq. (||) appears with positivity 



' " ~ 4 ) ~ 77 lo ^ ^' 



(70) 



8 



drj dV s pq ^/K 



r\ — x dr\ V x y/b. 



(71) 



Here soft band edge b n = a\ n . 

Equations obtained above are the starting point of the 
further analysis. 



B. Two-point kernel and universal eigenvalue 
correlations 



Two-point kernel determined by Eq. (23) can 
be calculated provided "wave function" ip n (x) — 
exp (—V s (x)) S n (x). Substitution of Eq. do9| ) into Eq. 
(E3) yields in the large-n limit 



K n (x,y) = 



IT b 



K n {x, y) = - kv ' 



(V-X) (xy) l/i {[b n -x] [6 n -y]} 1/4 



x {cos cf>„_i (x) cos $„ (y) - cos <I> n _i (y) cos <f>„ (x] 



(72) 



if x and y lie within the band (0, b n ). If at least one of 
the arguments in two-point kernel is negative, it is iden- 
tically zero (due to presence of hard edge) . In Eq. ( |7^ ) 
k n stands for leading coefficient of S n (x) . 
Taking into account the large-n identity 



&n-i (x) = $ n (x) — 2 arccos 
we obtain 

1 



* K (y- x )(xy) 1/4 {{b n ~x}[b n -y}} 1/4 



(73) 



x | cos (x) cos (y) - sin (y) cos $„ (x) \ ~ \ 



(74) 



Smoothed (over the rapid oscillations) connected corre- 
lator v c (x, y) of the density of eigenvalues, Eq. 

v c (x, y) = 



K (x + y)/2- xy 



27T ( X - y) y/xyy/bn - Xy/b~n~ 



!J 



x^y (75) 



manifests dependence on the potential V (x) only 
through the soft edge b n of the spectrum. 

The local properties of two-point kernel are obtained 
by assuming that in Eq. ( [74] ) \x — y\ <C b n and both x 
and y stay away from the hard edge x = and soft edge 
x = b n : 



K n {x, y) 



sin (irjvn bn (Qdg) 
7r (y - x) 



(76) 



The characteristic scale of the changing of f2b re (£) is of the 
order of b n , so that for \x — y\ <C b n Eq. ( f7(f ) is reduced 
to universal form Eq. (|^|) with V n = f2 bii (^j^) playing 
the role of local density of levels. Correspondingly, lo- 
cal two- level cluster function Y 2 (s, s') being rewritten in 
rescaled variables s and s' follows the universal form Eq. 
( |35j ) that proves universal Wigner-Dyson level statistics 
in the bulk of the spectrum for unitary random-matrix 
ensembles with confinement potentials V s (x). 



C. Density of levels and one-point Green's function 



Density of levels is obtained from Eq. ([76|) in the limit 

y -> x: 



K (x)) = \v 



drj dV s pq \/b n — x 



Ti" 2 Jo V ~ x drj V x y/b„ — r\ 



(77) 



Considering this expression as an equation for dV s /dx one 
can resolve it [ p6| arriving to the mean-field equation by 
Dyson 

V s (x)= [ dx' {v n (x'))ln\x-x'\+ fj,, (78) 



where integration runs over x' € (0, b n ). We once more 
stress that the mean-field equation is a direct conse- 
quence of the point-wise asymptotics for corresponding 
orthogonal polynomials S n {x) which involve Szego func- 
tion as a starting point. 

Correspondingly, the one-point Green's function 



G p (x) 



dV s _ f ™ drj dV s frj y/b n - x 

dx 7r J rj — x drj V x yjb n — rj 



(79) 
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D. Connected two-point Green's function 

In the maximum entropy models the smoothed con- 
nected two-point Green's function can be calculated in 
the same way as it was done in Section V. The only dif- 
ference is that integrals in Eqs. (EtI) - (|4^), ( [54] ) now run 
from to b n . Carrying out this integration with two- 
point kernel K n (x, y) given by Eq. (74) we arrive at the 
universal formula 



bn 
2 



(X- 



V 



(80) 



In contrast to the smoothed connected two-point 
Green's function the local one is determined by the same 
formulae Eqs. ( p0[ ) - (^2|) provided x and y are far from 
both edges. 



VIII. CONCLUSION 

We have presented rigorous analytical consideration 
of the matrix model given by non-Gaussian distribu- 
tion function P({x}), Eq. (Q), with very general class 
of confinement potentials V (x) within the framework 
of orthogonal-polynomials technique. Our treatment is 
equally applied to the random matrix models with pres- 
ence and absence of the hard edge in the eigenvalue 
spectrum. We have calculated with asymptotic accuracy 
the density of levels, one-point Green's function, two- 
point kernel, "density-density" correlator and two-point 
Green's function over all distance scale. 

It was established that two-point correlators in consid- 
ered random-matrix model posses a high degree of uni- 
versality. In the absence of hard edge the universality 



is observed for very wide class of monotonous confine- 
ment potentials V (x) which behave at least as 
(6 > 0) and can grow as or even faster than any poly- 
nomial at infinity (the case of border level confinement 
when V (x) ~ |a;| as \x\ — > oo has been treated in Ref. 
p2j). In the presence of hard edge in eigenvalue spec- 
trum the universality holds for monotonous confinement 

potentials V s (x) which behave at least as |a;| 5+<5 (S > 0) 
and can grow faster than any polynomial at infinity. 

We have shown that in those unitary non-Gaussian 
random-matrix models the density of levels and one- 
point Green's function essentially depend on the mea- 
sure, i. e. on the explicit form of confinement po- 
tential. In contrast, (connected) two-point characteris- 
tics of spectrum ("density-density" correlator, two-point 
Green's function) are rather universal. Indeed, we have 
observed global universality of smoothed two-point con- 
nected correlators and local universality of those with- 
out smoothing over rapid oscillations. In both cases the 
correlators were shown to depend on the measure only 
through the endpoints of spectrum (global universality) 
or through the local density of levels (local universality). 

Rigorous polynomial analysis enabled us to recover the 
results obtained before by different approximate methods 
and to extend previously known results for much wider 
class of random-matrix ensembles with strong confine- 
ment potentials irrespective of presence/absence of hard 
edge. We also have established a new local universal rela- 
tionship in the random- matrix theory for normalized and 
rescaled connected two-point Green's function gP p (s, s') 
[see Eq. (|62|)]. Finally, it is worthy of notice an interest- 
ing and quite surprising intimate connection between the 
structure of Szego function and mean-field equation that 
has been revealed in the proposed formalism. 
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Appendix: Integral representation of $„ (x) 

To prove Eq. (|3l]) let us calculate the first derivative of 7„ (x) (the calculations are similar to those done in Ref. 
!7|, Ch. 11). From Eq. ([H]) we obtain 

d ln _ i x _ v f +a - & (Ode 



dx 2n^T^ J_ an 

+ j^r - ^ „ (A1) 



whence 



,/ a 2 -^^iH = _L-p 



_ i /• +a " / e , 



2^ 7- a „ ^e i e-z 

After integration by parts we have 



^ 2tt y_ 0n £-x ^ - £ 2 £-3 y 

^TMO^f^V (A2) 



" 7n - 1 - r >#^ef- (as) 



^ tt^/o 2 - x 2 7 - * 2 d£ 

Substituting Eq. dn3) into Eq. (A3) and using identity 



p fTOFF = (A4) 



we obtain 



din 2 [ an ^ V<EIl t dV 1 / A r\ 



dx ^a 2 - x 2 Jq e 2 -a; 2 d£ 2^ 

The integral in Eq. (A5) may be handled as follows 



v r % = v r e« ^ v^i 



e 2 - x 2 * de y e 2 - a; 2 de v^e 1 
ede dv 



Bearing in mind Eq. ( jll| ) and introducing function 



2l v r & dv v° 

^ Jo e - a; 2 de 



the derivative d"f n /dx can be rewritten as 



dx 

Further, noting from Eq. (^ll) that j n (0) = 0, we obtain the integral representation 



a. 



or, equivalently (see Eq. §21 

$„ (x) = i arccos f - tt ^ cj q „ (£) d£ + ^ (2n - 1) . 



(A6) 



"«S*) = ^V\ ^ 2 ^^== (AT) 



\/a n — x 



7„ (x) = -7T / w Qri (e) d£ + In- - ) arcsin ( — ) , (A9) 



(A10) 
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